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ABSTRACT 

We summarize the studies and discussions of the Structure 
Function subgroup of the QCD working group of the Snowmass 
1996 Workshop; New Directions for High Energy Physics. 

I. INTRODUCTION 

Our knowledge of the structure functions of hadrons, and the 
parton density functions (PDFs) derived from them, has im¬ 
proved over time, due both to the steadily increasing quantity 
and precision of a wide variety of measurements, and a more 
sophisticated theoretical understanding of QCD. Structure func¬ 
tions, and PDFs, play a dual role; they are a necessary input 
to predictions for high momentum transfer processes involv¬ 
ing hadrons, and they contain important information themselves 
about the underlying physics of hadrons. Their study is an es¬ 
sential element for future progress in the understanding of fun¬ 
damental particles and interactions. 

Because of the ubiquitousness of structure functions, the ac¬ 
tivities of the subgroup had significant and productive overlap 
with several other subgroups, and were focused in a number of 
different directions. This summary roughly follows these di¬ 
rections. We start with the precision of our knowledge of the 
PDFs. There was an attempt to define a ‘Snowmass conven¬ 
tion’ on PDF errors, reviewing the experimental and theoretical 
input to the extraction of the PDFs and an appraisal of what is 
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left to do. Next, we explore the important connection between 
the strong coupling constant, ag, and the structure functions. 
One of the important inputs provided by the structure functions 
is in the precise extraction of electroweak parameters at hadron 
colliders. The systematic uncertainties in the structure func¬ 
tions may be the limiting factor in the determination of elec¬ 
troweak parameters, and this is discussed in the subsequent sec¬ 
tion. There is then a review of some relevant aspects of heavy 
quark hadroproduction. Finally, as a summary, we present an 
{x, Q^} map of what is known and what is to come. 

II. PRECISION OE PDES AND GEOBAE 
ANAEYSES 

The extraction of PDFs from measurements is a complex pro¬ 
cess, involving information from different experiments and a 
range of phenomenological and theoretical input. 

A. Experimental systematie errors 

Since the extraction of the PDFs usually requires using data 
from different experiments, and since the most precise data are 
usually limited by systematic, rather than statistical, errors, it 
is important that the systematic errors are taken properly into 
account. In particular, it is necessary to understand the correla¬ 
tions of different systematic errors on the measurements within 
and across experiments. Several groups have begun to make this 
information available in electronic and tabular form. Contribu¬ 
tions to these proceedings by Tim Bolton (NuTeV) and Allen 
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Figure 1: Fixed target DIS data. Note the full {x, Q^} region is 
clipped by the plot. 
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Figure 3; ep collider data. Note the full {cc, Q^} region is 
clipped by the plot. 
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Figure 2: Drell-Yan (E605), Direct Photon (E706, WA70, Hadron-hadron collider data. 

UA6), and DY asymmetry (NA51) data. 


Caldwell (ZEUS) give the details. 

B. {x, Q^} Kinematic Map for PDFs 

Global QCD analysis of lepton-hadron and hadron-hadron 
processes has made steady progress in testing the consistency 
of perturbative QCD (pQCD) within many different sets of data, 
and in yielding increasingly detailed information on the univer¬ 
sal parton distributions]^ We present a detailed compilation of 
the kinematic ranges covered by selected experiments from all 
high energy processes relevant for the determination of the uni¬ 
versal parton distributions. This allows an overall view of the 
overlaps and the gaps in the systematic determination of parton 
distributions; hence, this compilation provides a useful guide to 
the planning of future experiments and to the design of strate¬ 
gies for global analyses. 

These analyses incorporate diverse data sets including fixed- 


target deeply-inelastic scattering (DIS) data[^ |3]| of BCDMS, 
CCER, NMC, E665; collider DIS data of HI, ZEUS; lepton 
pair production (DY) data of E605, CDF; direct photon data 
of E706, WA70, UA6, CDF; DY asymmetry data of NA51; 
W-lepton asymmetry data of CDF; and hadronic inclusive jet 
datalQ] of CDF and DO. The total number of data points from 
these experiments is ^ 1300, and these cover a wide region in 
the kinematic {x, Q^} space anchored by HERA data at small 
X and Tevatron jet data at high Q. 

We now present the various experimental processes. Note that 
while this is a comprehensive selection of experiments, it is by 
no means exhaustive; we have attempted to display those data 
which are characteristic for the structure function determina¬ 
tion. In some cases, we have taken the liberty to interpret the 
data so as to facilitate comparison among the diverse processes 
we consider.^ Also note that we have not attempted to deal with 
the different precision of different measurements, or to sepa¬ 
rately consider the quark and gluon determination; the reader 


^PDF sets are available via WWW on the CTEQ page at 


htBi://www.phys.psu.edu/'-s.cteqj and on the The Durham/RAL HEP Database 
at |http://durpdg.dur.ac.uk/HEPE)ATA/HEPDATA.htm] . 


^In particular, since we have taken the data points from the global fitting 
files, there is a cut on the minimum value of Q 2GeV to avoid the non- 
perturbative region. 
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Figure 5; Comparison of gluons obtained with pre-1995 DIS 


data (A-series) with those using current DIS data (B-series). 
(Cf., Ref. %.) 
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Figure 6: Comparison of gluons obtained without jet data (C- 


series) with those obtained with jet data, D-series (CTEQ4Ax). 
(Cf., Ref. [0.) 


should keep these points in mind when comparing the figures. 

The quark distributions inside the nucleon have been quite 
well determined from precise DIS and other processes, cf., 

Fig.0 

+ N ^ + X . (1) 

Improved DIS data in the small-x region is available from 
FIERA, and this is of sufficient precision to be sensitive to the 
indirect influence of gluons via high order processes. 

The Drell-Yan process is related by crossing to DIS. In lowest 
order QCD it is described by quark-antiquark annihilation; 

q + q^-f* ^ e++r . ( 2 ) 

The kinematic coverage is shown in Eig. ^ and Eig. ^ 

Recent emphasis has focused on the more elusive gluon distri¬ 
bution, G{x, Q), which is strongly coupled to the measurement 
of as- Direct photon production, 

g + q^i + q , q + q^i + g , ( 3 ) 

in particular from the high statistics fixed target experiments, 
has long been regarded as the most useful source of information 
on G{x, Q), cf., Eig. However, there are a number of large 
theoretical uncertainties (e.g., significant scale dependence, and 


kx broadening of initial state partons due to gluon radiation)!^ 
^ that need to be brought under control before direct photon 
data can place a tight constraint on the gluon distribution. 

Eor example, the broadening due to soft gluon radiation 
is essentially a higher twist effect (but with a large coefficient), 
and should affect all hard scattering cross sections. The mag¬ 
nitude of the correction to the cross section should be on the 
order of n{n + l)(fc'r)^/(4p|’), where (kr) is the average in 
the hard scatter, and n is the exponent of the differential cross 
section with respect to pT ■ (da/dpT oc XjpVf). Eor the Teva- 
tron collider regime, the effect should fall off as ~ 1 /p^, as 
is observed for example in direct photon production in CDF 
For pt > 50 GeV, the effect is negligible. For fixed target ex¬ 
periments, the effective value of n is large and changes rapidly 
with Pt (due to the rapidly falling parton distributions). The 
soft gluon radiation tends to make the cross section steeper at 
low Pt and at high pT, and to cause an overall normalization 
shift of a factor of 2.[^ There are several approximate methods 
to predict the effects of soft gluon radiation, as for example in 
gaussian smearing, or the incorporation of parton showers 
into a NLO Monte Carlo. Further understanding may await the 
development of a more formal treatment of the effect. Several 
theoretical ideas are under development. 

Inclusive jet production in hadron-hadron collisions, 

{gg, qq} {gg, qq} , g + {q, q} ^ g + { 9 , q}, (4) 

is very sensitive to as and G{x, Q), (Fig. NLO inclusive 
jet cross sections yield relatively small p scale dependence for 
moderate to large Et values. 0 High precision data on single 
jet production is now available over a wide range of energies, 
15 GeV < Et< 450 GeF.0 Eor Et > ^QGeV, both the the¬ 
oretical and experimental systematic errors are felt to be under 
control. Thus, it is natural to incorporate inclusive jet data in a 
global QCD analysis. 

In reviewing the figures we see the large kinematic range 
which is explored by these processes. It is a useful exercise 
to overlay the curves according to the separate determination 
of the valence-quarks, light-sea-quarks, heavy-quarks, and glu¬ 
ons. Although there is no room here for such a presentation, 
we leave this as an exercise to the interested reader. Obviously, 
when comparing such a wide range of processes, one must keep 
in mind considerations beyond just the kinematic ranges. Eor 
example, the DIS and Drell-Yan processes are useful in deter¬ 
mining the quark distributions, whereas the direct photon and 
photoproduction experiments yield information about the gluon 
distributions-though not with comparable accuracy; the deter¬ 
mination of the gluon distribution is subject to many more the¬ 
oretical and experimental uncertainties. 

Likewise, the systematics for hadron-hadron and lepton- 
hadron processes are quite different. Specifically, while the 
hadron-hadron colliders can in principle determine parton dis¬ 
tributions out to large Q^, extractions of PDEs from these data 
are only beginning. DIS experiments probe small x (HERA) 
and high x (NuTeV), and low-mass Drell-Yan collider measure¬ 
ments yield complementary results at higher Q^. This combina¬ 
tion of experiments improves the reliability of the PDEs, allows 
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for cross checks among the different experiments, and yields 
precise tests of the QCD evolution of the parton distributions. 

1. Progress of PDFs 

As new global PDF hts are being updated and improved, it 
can be difficult to quantify our progress as to how precisely we 
are measuring the hadronic structure. To illustrate this progress 
we consider sets of global PDF hts using various subsets of the 
complete data set.^ 

First, we compare the A- and B-series of hts shown in Fig. ||. 
The A-series shows a selection of gluon PDFs extracted from 
pre-1995 DIS data using various values of as as indicated 
in the hgure. The B-series shows the same selection, but in¬ 
cluding the recent DIS data. By comparing the A- and B-series 
of hts, we found that recent DIS data of NMC, E665, HI 
and ZEUS considerably narrow down the allowed range of the 
parton distributions. 

Next, we compare the B- and C-series of hts shown in Eig. |] 
and Eig. These hts were performed with the same data set, but 
the C-series ht used a more generalized parametrization with ad¬ 
ditional degrees of freedom. By contrasting the B- and C-series 
we see that we must be careful to ensure that our parameteri¬ 
zation of the initial PDEs at Qo is not restricting the extracted 
distributions. 

Einally, we compare the C- and D-series (CTEQ4Ax) of hts 
shown in Eig. |[ Eor the D-series hts, the Tevatron jet data was 
used, whereas this was excluded from the C-series hts. The jet 
data has a signihcant effect in more fully constraining G{x, Q) 
as compared to the C-series. The quality of the hnal D-series 
hts (CTEQ4Ax) is indicative of the progress that has been made 
in this latest round of global analysis. 

C. High Ei Jets and Parton Distributions 




Eigure 7: CDE and DO data compared to NLO QCD using a) 
CTEQ4M and b) CTEQ4HJ. C/, Refs. 


High statistics inclusive jet production measurements at the 
Tevatron have received much attention recently because the high 
jet Et cross-section[p^ [T^ Q is larger than expected from 
NLO QCD calculations. A comparison of the inclusive jet 
data of CDE and DO and results is given in Eig. We see that 

^For the details of how these fits were performed, see the original paper. 
Ref. IQ]. 



Gluon X 


Eigure 8; (a) The CTEQ4HJ gluon distributions are compared 
to that of CTEQ3M: (b) the ratio of the CTEQ4HJ gluons to 
CTEQ3M. Cf, Ref. [g]. 


Table I: Total /point) values and their distribution among 

the DIS and DY experiments for CTEQ4M and CTEQ4HJ. Cf, 

Ref. [ 0 ]. 


Experiment 

#pts 

CTEQ4M 

CTEQ4HJ 

DIS-Eixed Target 

817 

855.2(1.05) 

884.3(1.08) 

DIS-HERA 

351 

362.3(1.03) 

352.9(1.01) 

DY rel. 

129 

102.6(0.80) 

105.5(0.82) 

Total 

1297 

1320 

1343 
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there is a discernible rise of the data above the fit curve (hor¬ 
izontal axis) in the high Et region. The essential question is 
whether the high Et jet data can be explained in the conven¬ 
tional theoretical framework, or require the presence of “new 
physics”. [010 

Although inclusive jet data was included in the global fit 
of the PDF, it is understandable why the new parton distri¬ 
butions (e.g., CTEQ4M) still underestimate the experimental 
cross-section: these data points have large errors, so they do 
not carry much statistical weight in the fitting process, and the 
simple (unsigned) total is not sensitive to the pattern that the 
points are uniformly higher in the large Et region. A recent 
study investigated the feasibility of accommodating these data 
in the conventional QCD framework by exploiting the flexibility 
of G{x, Q) at higher values of x (where there are few indepen¬ 
dent constraints), while maintaining good agreement with other 
data sets in the global analysis. [Q] 

A result of this study is the CTEQ4HJ parton sets which are 
tailored to accommodate the high Et (> 200 GeV) jets,^ as 
well as the other data of the global flt.[Q] Eig. ^ compares pre¬ 
dictions of CTEQ4HJ with the results of both CDE and DO.^ 
Results shown in Eig. ^ and Table | quantifies the changes in 
values due to the requirement of fitting the high Et jets. Com¬ 
pared to the best fit CTEQ4M, the overall for CTEQ4HJ is 
indeed only slightly higher. Thus the price for accommo¬ 
dating the high Et jets is negligible. 

The much discussed high Et inclusive jet cross-section has 
been shown to be compatible with all existing data within the 
framework of conventional pQCD provided flexibility is given 
to the non-perturbative gluon distribution shape in the large-x 
region. 

Presently, we note that the direct photon data from the Eer- 
milab experiment E706 are sensitive to the same x range that 
affect the Tevatron high Et jet data. A more quantitative theo¬ 
retical treatment of soft gluons may allow the direct photon data 
to probe this question more precisely. We will need such accu¬ 
rate, independent measurements of the large-a; gluons to verify 
if the high-Et jet puzzle is resolved, or whether we have only 
absorbed the “new physics” into the PDEs. 

Nevertheless, this episode provides an instructive lesson: the 
precision with which we know the PDEs is not indicated from 
a simple comparison of different global fit sets. These fits pro¬ 
ceed from similar assumptions and procedures, so their relative 
agreement should not be taken as assurance of our knowledge 
of the PDEs. In the present case, the gluon density was naively 
estimated to be less than 10-20% (in the x kinematic range rel¬ 
evant for high Et jet production) from a simple comparison of 
different PDE sets. Surprisingly, a large change was eventually 
required (and accommodated) by the data (assuming the Teva¬ 
tron result is not an indication of some new physics). 


^This set is tailored to accomrnpdate the high Et jets by artificially decreas¬ 
ing the errors in the fit. See Ref. |M] for details. The of Table I is computed 
using the true en'ors. 

^For this comparison, an overall normalization factor of 1.01(0.98) for the 
CDF(DO) data set is found to be optimal in bringing agreement between theory 
and experiment. 


D. Challenges for Global Fitting 

Global fitting of PDEs is a highly complex procedure which is 
both an art and a science. This requires fitting a large number of 
data points from diverse experiments with differing systematics. 
Eurthermore, the data are compared with NLO theory which 
introduces additional complications on the theoretical side. 

There was extensive discussion as to how to determine the 
uncertainty of the PDEs. We note that one of the most impor¬ 
tant uncertainties for the PDEs is the choice of as, since this af¬ 
fects the gluon distribution directly as well as the singlet quarks. 
Both MRS and CTEQ now provide different PDEs with differ¬ 
ent choices of as, this is a big improvement in determining the 
uncertainty of PDEs. But this group did not succeed in answer¬ 
ing all the questions related to the goal of a true one-standard- 
deviation covariance matrix of PDE uncertainties, although we 
did focus on some points that deserve further study. We list 
some of these below. 

1. A reminder: when adding two experiments, you simply 
add their x^’s, and Ax^ = 1 of the total is one a. 

(a) Due to direct photon theory /r-scale and kr uncer¬ 
tainties, there is no way to define one standard devi¬ 
ation for these data. The handling of the /r-scale is 
done differently in different groups and can lead to 
somewhat different gluon distributions. 

(b) Other “choices” can lead to significant differences in 
X^ { Ay^ « 50 — 100 units is typical). These choices 
include which data sets to use, the starting Qq value, 
etc. One example is the small-x CCER neutrino data 
which disagrees with the electron/muon DIS data. 
This difference is unlikely to be caused entirely by 
parton distributions, and how this is handled in the 
global fits can cause significant changes in the global 

2. Many experiments do not provide correlation matrices, and 
we’ve never seen a correlation matrix for a theoretical un¬ 
certainty. Without both of these for every experiment, one 
cannot expect Ax^ = 1 to work. 

3. In principle we should add in LEP/tau/lattice constraints on 
as- But if they are treated as only a single data point, they 
will be swamped by the other 1200 points. (This would not 
be true if Ax^ = 1 were valid.) 

4. What to do about the charm mass in DIS? It will change E 2 
predictions, but the resulting parton distributions then have 
a different definition of “heavy quark in the proton” and 
this must be accounted for in the theoretical calculations. 

5. When the CDE W asymmetry and NA51 data were added 
(the change between CTEQ2 and CTEQ3), they gave a 
consistent picture of u and d. But the went up for the 
rest of the experiments by 30! Once again Ax^ = 1 is 
invalid. The choice was to accept the larger x^’s to incor¬ 
porate the new data presented. 
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6. What about higher twists? Should a higher twist theoreti¬ 
cal uncertainty be added to DIS data? 


E. Choice of Parametrization 




Figure 12; The ratio of a) c{x)ls{x) and b) h{x)/c(x) for a 
range of Q. For increasing Q, the evolution reduces any differ¬ 
ence between the distributions. 


Figure 9: a) The gluon distribution xG{x) at Qq = 1.6GeV 
using the MRS and CTEQ parametrizations. The two curves are 
indistinguishable in this plot, b) Fractional deviation for gluon 
of the CTEQ and MRS parametrizations. Note the full range of 
the y-axis is ±1%. 



Eigure 10; a) The u-valence distribution xuv(x) at Qq = 
l.QGeV using the MRS and CTEQ parametrizations. The two 
curves are indistinguishable in this plot, b) Eractional deviation 
for u-valence of the CTEQ and MRS parametrizations. Note the 
full range of the y-axis is ±1%. 



Eigure 11; a) The d-valence distribution xdy{x) at Qq = 
l.&GeV using the MRS and CTEQ parametrizations. The two 
curves are indistinguishable in this plot, b) Eractional deviation 
for d-valence of the CTEQ and MRS parametrizations. Note the 
full range of the y-axis is ±2%. 


The choice of boundary conditions for the PDEs at the ini¬ 
tial Qq has received increasing attention as the accuracy of data 
improves. [j^ ^ ^ Although the DGLAP evolution equation 
clearly tells us how to relate PDEs at differing scales, the form 
of the distribution at Qq cannot yet be derived from first prin¬ 
ciples, and must be extracted from data. Eor this purpose, it is 
practical to choose a parametrization for the PDEs at the initial 
Qq with a small number of free parameters that can be fit to the 
data. 


A question that was repeatedly raised in the workshop is the 
extent to which the choice of parametrization limits the ex¬ 
tracted PDEs of the global fit. It is important to note that the 
evolution equation for the global fits is solved numerically on 
an {x, Q^} grid. Therefore the issue of the parametrization is 
only relevant at Qq. Eor Q > Qq, the parametrized form is 
replaced by a discrete {x, Q^} grid. 

To approach these questions in a quantitative manner, we per¬ 
formed a simple exercise to examine the potential difference of 
PDEs that can be attributed to different choices of parametriza¬ 
tions. Specifically, we investigated the difference between the 
MRS0 and CTEQ0 parametrizations, which take the general 
form; 

MRS; 


f{x, Q) = aQX°‘^ (1 — (1 -f aQ^/x + a 4 x) 


(5) 


CTEQ; 


f{x,Q) = hQX^^{l - -f h^x'’*) 


^2 / 




(6) 


We used the CTEQ3M PDE set at Qq = 1.Q GeV (which is nat¬ 
urally described by the CTEQ parametrization shown above), 
and performed a fit to describe the same PDE set using the MRS 
parametrization. Note that this is an academic exercise that does 
not fit data, but rather explores the flexibility of the parametriza¬ 
tions. 

If we can accurately describe the CTEQ3 PDEs with the 
MRS form, then it is plausible that the particular parametriza¬ 
tion choice for the PDEs has little consequence. However, if 
we cannot accurately describe the CTEQ3 PDEs with the MRS 
form, we will need to investigate more thoroughly whether the 
parametrization introduces a strong bias as to the possible PDEs 
which will come from the global fitting procedure. 

In Eigs. H 1^ and 11, we plot both the CTEQ3 PDEs and 
the fit to the CTEQ3 PDEs using the MRS parametrized form, 
Eq.| We only show the gluon, u-valence, and d-valence; the 
results for the sea-quarks will be similar to the gluon. Eirst we 
plot X f{x) at Qq = 1.6 GeV on a Log-Log scale. The two sep¬ 
arate curves are indistinguishable in this plot. To better illus¬ 
trate the differences, we plot the fractional difference between 
the two PDEs. Observing that the scale on this plot is < 2%, we 
see that the variation over the range x = 1] is relatively 

small. We And larger deviations in the high x region, but the 
significance of this is diminished by the fact that the PDEs are 
small in this region. 
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Although we do not claim that this is an exhaustive investi¬ 
gation, this simple exercise appears to indicate that the PDFs 
extracted from a global fit should be insensitive to the choice 
of the above parametrizations (Eq. |]). One might speculate that 
the same conclusion would also hold for other parametrizations; 
however, such an exercise has yet to be performed. 

Furthermore, since the QCD evolution is stable as one evolves 
up to higher values of Q, any small differences at Qq will de¬ 
crease for Q > Qq. We can roughly see this effect by examing 
the ratios of c{x)/s{x) and b(x)/c[x) as shown in Fig. 0 For 
example, at Q = 10 GeV, the b-quark is less than half of the 
c-quark distribution; at Q = 100 GeV the b-quark distribution 
is significantly closer to the c-quark. [jll[] This observation sug¬ 
gests that the small differences we observed at Qq = 1.6 GeV 
will quickly wash out as we evolve upwards. 

The above observations, however, only apply in regions of 
X where PDFs are well-determined; and they cannot be taken 
literally without qualification. An important example which il¬ 
lustrates the importance of exercising caution is the behavior of 
the gluon distribution at large x brought to focus by the high 
Et jet data, as discussed in the last section. Whereas, using 
“conventional” parametrization of the form Eq.(|]), GMRS 0 
found it impossible to fit the jet results along with the rest of 
the global data, CTEQ showed that allowing for a more flexi¬ 
ble parametrization of the gluon distribution at large x can ac¬ 
commodate both. To accomplish this, one will need a func¬ 
tional form such as Eq.(||), with 64 substantially bigger than 1, 
or equivalently, in place of x^^. (Since 0 < x < 1, and 
the whole expression is multiplied by (1 — x)^^ which is steeply 
falling, G{x, Q) is still well-behaved.) The difference in the 
size of G{x, Q) resulting from these parametrizations can be as 
much as 100% at x = 0.5, as shown in Fig. |^. 

III. STRUCTURE FUNCTIONS AND as 

Structure functions are important in that they give us infor¬ 
mation on the value of as, and also in that they are often inputs 
to many different measurements, some of which themselves are 
used to determine a^,. In this summary of the work of the joint 
Og-structure function groups we investigate how structure func¬ 
tions themselves give us direct information on a^, and the ex¬ 
pected uncertainties of possible new measurements of structure 
functions at future proposed machines. There are two categories 
of structure function analyses which result in an Og measure¬ 
ment; Q^ evolution of structure functions, and measurements 
of sum rules, which pertain to the integrals of specihc structure 
functions over x. Since the theoretical and experimental errors 
are comparable for some of these analyses, this report examines 
how improvements might be made in both areas. 

On the experimental side of the study, we consider a /rp col¬ 
lider or an ep collider, and also a neutrino scattering experiment 
at a /r■*■/!“ collider. Given the current level of error in as mea¬ 
surements, we consider here only analyses which may result in 
few per cent or less error on ag(Mf). To address the theoretical 
issues within the scope of this report we can at best point out the 
largest problems and how they are currently being investigated, 
in hopes of inspiring theorists to devote more time to them. 


A. Evolution of Structure Functions 

When looking at the Q^ evolution of structure functions, one 
can use the DGF-Altarelli-Parisi Equations to hnd ctg [|^. In 
the case of the non-singlet structure functions the evolution as a 
function of Q^ is simply related to as and the non-singlet struc¬ 
ture function itself. In the case of the singlet structure functions, 
the evolution is related to ag, the structure function itself, AND 
the gluon distribution, which complicates matters. In either case 
care must be taken avoid large higher twist effects, which are 
particularly important at low Q^. 

Non-singlet structure functions can be measured in both neu¬ 
trino and charged lepton scattering experiments. One way is 
by taking the average of xF^^ and xF^^, where vN in the 
superscript indicates the presence of an isoscalar target. Simi¬ 
larly, averaging xW^'^ and xW^ also results in a pure non¬ 
singlet structure function, where the lepton is either an electron 
or muon, and the scattering center is a deuteron. Finally, one 
can use the structure function or F 2 ^ at high x, since 
there are virtually no sea quarks at high x. 

Many high statistics determinations of as have to date been 
performed, using a variety of techniques. By htting only xFq 
or F 2 at high x, one can do a pure non-singlet ht to the evo¬ 
lution, with no dependences on the gluon distribution. Given 
the wealth of precise data in charged lepton scattering structure 
functions, however, there are also determinations of as from 
htting F 2 at all x, but including a contribution to the evolution 
from the gluons. These two different kinds of determinations 
do not show aiw systematic difference in the hnal result, as is 
shown in table The errors listed in ^ are deceiving, however. 


Table II: A selection of as measurements from structure func¬ 
tions, and the total error on as{M‘^). 


Method Experiment 


xFq only 
xFq and F 2 
F 2 low X 
F 2 high X 
F 2 low X 


CCFR[D| 
CCFR [0 
NMC @ 


SFAC/BCDMS [jl^l 
HERA [0 


as{Ml) 


25 

.118 

± 

.007 

25 

.119 

± 

.0055 

7 

.118 

± 

.015 

50 

.113 

± 

.005 

4-100 

.120 

± 

.010 


because in fact they are all dominated by either experimental 
or theoretical systematic errors. In the remainder of this sec¬ 
tion we consider the largest two systematic errors, and how new 
machines (and new calculations) could hopefully reduce these 
errors. 

1. Experimental Errors on as and possible improvements 

The dominant experimental systematic error in the measure¬ 
ments listed in the table come from energy uncertainties. These 
can come from spectrometer resolution, calibration uncertainty 
in the detector, or overall detector energy scale. The key to im¬ 
proving the overall experimental error in these measurements is 
not higher statistics or higher energies, but better calorimetry, 
and better calibration techniques. The challenge in determin- 
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ing the energy scale in deep inelastic scattering experiments is 
in finding some “standard candle” from which to calibrate. For 
example, if there were some way of measuring the known mass 
of some particle decaying in the system, or if the initial beam 
energy was very well known because of accelerator constraints, 
this could substantially improve the energy scale determination 
over previous experiments. 

A number of machines have been proposed at this workshop 
in a variety of energies and initial particles. While it is true that 
machines (and experiments) are not proposed these days to do 
precise QCD measurements alone, there are some interesting 
possibilities that may arise from these machines. 

Because of other considerations (namely the rise of F 2 at low 
x) a lepton/hadron collider is an attractive possibility. Currently, 
however, the HERA as experimental error is dominated by un¬ 
certainties in the x distribution of the structure functions mea¬ 
sured, (particularly that of the gluon). In order to do a DGLAP- 
style evolution measurement in a lepton hadron collider, one 
would need to have both t^p and i~p collisions, measure the 
different cross sections, and extract a non-singlet structure func¬ 
tion. The statistics needed for a precise structure function mea¬ 
surement at the energies that have been proposed would be well 
above current HERA expectations, and the higher in these 
machines operate the lower the cross section, and the smaller 
the effect one is trying to measure. 

Another intriguing possibility would be a neutrino experiment 
at a muon collider. A 2TeV muon collider could (with consid¬ 
erable engineering) make very high rate SOOGeV neutrino and 
antineutrino beams. If one knew the muon beam energy very 
well (taking as an example how well the LEP energy scale is 
now known after much work!) then a neutrino beam coming 
from muon decays would be at a very well-understood energy 
as well. There would be negligible production uncertainty from 
a neutrino beam coming from a muon beam, and the rates for 
such a beam would be astronomical simply starting with the 
current proposals for muon intensities in the accelerator. 

2. Theoretical Errors on as and possible improvements 

Currently the renormalization and factorization scale uncer¬ 
tainties dominate the theoretical error on as from structure 
function evolution. This is true for both singlet and non-singlet 
structure functions evolution. By assuming the factorization and 
renormalization scales were and k 2 Q^ respectively, and 

varying ki and ^2 between 0.10 and 4, Virchaux and Milsz- 
tajn atTive at an error of 5{as{M‘^)) = .004. [p^. They claim 
that the overall of the fit did not increase signihcantly when 
these variations were made. Similar or larger QCD scale errors 
apply to the other as measurements listed in table Certainly 
the most straightforward (and perhaps naive) way to reduce this 
error would be to calculate the next higher-order term in the 
DGLAP equations. 

Still another method of reducing these errors is to actually ht 
for ki and /c 2 , and see what the resulting error on these values 
are within the fit. By floating those constants, however, one is 
assuming QCD works, and getting a good fit for one consistent 
value of A-pg in the experiment can no longer be claimed by 
itself as a test of QCD. If ki and k 2 are floated, one does not 


test QCD until one compares one experiment’s as value with 
another experiment’s value. Eurthermore, if the fit prefers val¬ 
ues of fci and k 2 far away from one, one would also question 
the validity of the measurement. 

B. Sum Rules are Better than Others 

The two sum rules that have thus far been used to measure as 
are the Gross Llewellyn Smith sum rule plj ] and the Bjorken 
sum rule [^, which are related to xF^, and the polarized struc¬ 
ture functions Pn(x) and gp{x) respectively. Since these meth¬ 
ods of determining as are far less mature than the structure 
function analysis, the corresponding experimental errors on as 
are much larger. Since both Sum Rules come are fundamental 
theoretical predictions, and the higher order corrections to the 
sum rules are so straightforward to compute, the QCD scale er¬ 
ror on these measurements is much smaller than those of the 
evolution measurements. Table n gives a list of systematic and 
statistical errors for both sum rules. 

Table 111: Table of errors on as{M‘^) from Sum Rules. ^ Erom 
E142 result, E154 claims should be higher 


Error Source 

5a 

GLS 

s{Ml) 

Bjorken 

Statistical 

.004 

< .001 

Low X extrapolation 

.002 


Overall Normalization 

.003 

.002 

Experimental Systematics 

.004 

.006 

Higher Twist 

.005 

.003-.008 

QCD Scale Dependence 

.001 

.002 


1. Low x Uncertainties 

The largest uncertainties in sum rule measurements come 
from the fact that they involve integrals from a; = 0 to x = 1. Of 
course no experiment can measure all the way down to x = 0, 
and the closer to 0 one can reach the smaller the error in ex¬ 
trapolating from the lowest data point to zero will be. What is 
usually done to extrapolate to x = 0 is a functional form is as¬ 
sumed, and the data are either ht to that functional form and the 
resulting parameters are checked with a theory, or if the data 
do not have enough statistical precision some functional form 
is simply assumed. While for the GLS sum rule the data seem 
to agree with simple quark counting arguments for the form of 
xFa at low x, the newest data from SLAC El54 (shown after 
Snowmass’96 at 1CHEP96) do not ht to a function whose inte¬ 
gral converges as x goes to 0. The collaboration does not yet 
report a measurement of as from their new data, saying that the 
low X behavior of the integral is too uncertain; however, this 
analysis is in progress. For future improvement on the Bjorken 
Sum Rule one will need to go to lower x than what has cur¬ 
rently been reached (x = .015). If the Bjorken integral is not 
hnite then much more is called into question than the validity of 
QCD! 
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Another uncertainty associated with the low x region is that 
one also needs to go to low to measure low x. At low 
higher twist uncertainties become important, and these higher 
twist contributions have never been measured for these sum 
rules. The present state of higher twist calculations for DIS 
sum rules is given in reference which discusses results 
from many models of higher twist calculations, including QCD 
Sum Rules, and a non-relativistic quark model. Again, there 
is more trouble associated with the Bjorken Sum Rule than the 
GLS sum rule, because the different models predict very dif¬ 
ferent higher twist contributions to the former, while agreeing 
at the 50% level for the latter. So, whether one takes as the 
higher twist error the spread of theoretical predictions or the er¬ 
ror on one such prediction (shown in the table above) one can 
arrive at very different errors. In either case that error is sig- 
nihcant at the cutTently relevant region. Unless a proven 
agreed-upon method of higher twist calculations arises the best 
bet in the future will be to simply ht the sum rule results for a 
higher twist contribution and an as contribution. This will re¬ 
quire much higher statistical precision in the structure function 
measurements themselves than what is cutTently available. 

2. Normalization Uncertainties 

Finally, if one proposes to improve these measurements by 
going to a higher the next most important error (assuming 
one has solved the problem of extrapolating to low x) will be 
the overall normalization error. Since the effect one is measur¬ 
ing is proportional to 1-ergand not as, as gets larger and as 
gets smaller then an overall 1% error on the normalization of 
the structure functions (and hence the integral itself) turns into 
a larger fractional error on as- This is shown quantitatively in 
hgure 1^ which shows the effects of the higher twist error as 
a function of and a 1% normalization error on the struc¬ 
ture functions as a function of Q^. The sum of the two etTors 
in quadrature show that measuring the sum rules at a above 
100 GeV^ will not reduce the overall error for even an ambitious 
normalization error of 0.25%. The current normalization error 
on the overall j/-nucleon cross section is 1%, and the error on 
the ratio of v and V cross sections is another 1%, which trans¬ 
lates into presently a total xF^ normalization of 1.4%. There 
are currently no plans to improve this measurement, one would 
need a tagged neutrino beam (which might be possible with a 
muon beam at a muon collider) to do so. 

C. as conclusions 

Structure functions and QCD provide us with the possibility 
of two complementary measurements of as', the evolution 
and sum rules. The curi'ent etTors on as from structure function 
evolution are in the 4 — 5% range, and will be improved only 
with a reduction of the renormalization and factorization scale 
uncertainties. For this, next to next to leading order (NNLO) 
corrections to the DGLAP equations must be computed. By far 
the most important experimental uncertainty in evolution mea¬ 
surements comes from how precisely experiments know their 
energy scale and resolution. Sum Rules have very different out¬ 
standing issues, namely the low x uncertainty and measurement. 


and also the higher twist terms. The best way to eliminate higher 
twist uncertainties would be to simply measure their contribu¬ 
tions in the lowest yet still have enough statistics at higher 
for a measurement of as- While the sum rule analyses would 
beneht from much higher statistics, in general, to arrive at new 
measurements of as from structure functions we must do more 
than simply raise the energies of the experiments and run them 
longer! 

IV. STRUCTURE FUNCTION INPUTS TO 
PRECISION EEECTROWEAK 
MEASUREMENTS 

Structure functions are inputs to many precision electroweak 
measurements-a few examples are from sin^Ow measured in 
r/N scattering and global electroweak hts which include as from 
structure function data along with other fundamental parame¬ 
ters. A measurement expected to have signiheant experimental 
improvement in the future such that the structure function uncer¬ 
tainty becomes important relative to other uncertainties is the W 
mass (Mw) measurement from on-shell production at collider 
experiments. Even at the current level of precision of this anal¬ 
ysis there are outstanding questions about how that uncertainty 
is evaluated, and whether this could be improved, even before 
new experiments come around. 

At a hadron collider experiment, the W mass itself cannot 
be directly measured on an event by event basis, because the 
clean signatures of W production contain a charged lepton 
and therefore also contain a neutrino. Furthermore, the initial 
center of mass energy of the partons which interact to give a 
W is not known, so one cannot simply require the total mo¬ 
menta to balance to give the energy of the outgoing neutrino. 
One can use the constraint that the total initial transverse mo¬ 
mentum is zero, however. The way the mass is then mea¬ 
sured in an experiment is that the transverse mass is computed, 
Mt = y/2pk(l — cos where are the transverse mo- 
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Figure 13: Variation of (a) higher twist, (b) 1% normalization 
errors, and (c) the sum in quadrature of the two as a function of 
for the GLS sum rule. 
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menta of the charged lepton and neutrino, and (j)^" is the angle 
between the charged lepton and neutrino in the transverse plane. 
The shape of the Mt distribution is then extremely sensitive to 
Mw, but is also dependent on the parton distributions used in 
the Monte Carlo simulation, in particular, the transverse com¬ 
ponent of u — d. 

Table gives the uncertainty in Mw from CDF and DO from 
direct production, and measurement of the transverse mass . 

Currently the structure function uncertainty is estimated by 
doing the analysis with several different sets of parton distribu¬ 
tion functions, and comparing the results, using the W asym¬ 
metry measurement as another constraint. Figure shows the 
measured W asymmetry from CDF and the predictions from 
various PDFs Given that most of these PDFs come from 
the same input data (deep inelastic structure functions), the 
spread of the predictions represents an error in the technique of 
parametrizing the distribution which accounts for the W asym¬ 
metry, not the error on the distribution itself. By requiring a 
PDF to reproduce the measured W asymmetry, one is choosing 
a more appropriate parametrization, but one must go further to 
assign errors on that specihc parametrization. 

Figure shows the resulting change in Mw for different 
PDFs, and how many standard deviations each PDF is from pre¬ 
dicting the W asymmetry [ psj . 

The problem with estimating this uncertainty by comparing 
different PDFs is the following: if all of these PDFs are sim¬ 
ply different parametrizations which come from the same sets 
of deep inelastic scattering data, then two different PDFs do not 
necessarily encompass the uncertainty on whatever quark distri¬ 
butions are relevant. There must be errors on the PDFs in order 
for the correct error on the W mass uncertainty to be evaluated. 


Table IV: Table of uncertainties for both the CDF and DO W 
mass measurements, for different hnal states (ei/ or piz) and 
different run periods (Ia,Ib). 


Source 

CDE 

DO 

DO 



la 

la 

Ib 


e 


e 

e 

Statistics 

145 

205 

140 

70 

Lepton Scale 

120 

50 

160 

80 

Lepton Resolution 

80 

60 

85 

50 

Lepton Efficiency 

25 

10 

30 

20 

Py, PDF 

65 

65 

65 

65 

pKecoil ]y[o(jgj 

60 

60 

100 

55 

Underlying Event 
in Lepton Towers 

10 

5 

35 

30 

Background 

10 

25 

35 

15 

Trigger Bias 

0 

25 

- 

- 

QCD Higher Order Terms 

20 

20 

- 

- 

QED Radiative Corrections 

20 

20 

20 

20 

Luminosity Dependence 

- 

- 

- 

70 

Total 


180 

270 

170 



0.0 0.5 1.0 1.5 2.0 

|Lepton Rapidityl 


Figure 14: The W asymmetry as measured in CDF and the pre¬ 
diction of various different parton distribution functions. 




Figure 15: The change in W mass versus the signed standard 
deviation of agreement with the measured W charge asymmetry 
for different PDFs. 


Of course, since at the present time there are no errors given 
with PDFs, this is not possible. 

There was much discussion at Snowmass about the difficul¬ 
ties associated with assigning errors to PDFs, and one should 
refer to that section of this write-up, and a separate submission 
by Tim Bolton on this topic. Given that the job of assigning 
those errors is one that is far from completion, a temporary so¬ 
lution was suggested at this meeting. Namely, a PDF-generator 
could produce a set of PDFs that span the range of the possible 
values of the distribution in question. For example, for the jet 
Et analysis, a set of PDFs with different values of as has been 
generated. Similarly, a set of PDFs with the acceptable range of 
u—d, which is important for the W mass measurement ( and also 
the W asymmetry measurement) could also be generated. Then 
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the Mw analysis could simply compare the different PDFs in 
one set provided for an estimate on the Mw error from uncer¬ 
tainty in the PDFs, and compare different parametrizations for 
the uncertainty on the parametrization. 

Furthermore, much care must taken when using the measured 
W asymmetry to constrain the W mass error from PDFs. Since 
the W mass is measured using a distribution dependent mostly 
on transverse quantities, and the asymmetries depend on lon¬ 
gitudinal differences between the u and d quark distributions, 
the correlations (and/or lack of correlations) must be taken into 
account appropriately. 

In order to use PDFs to their full potential, and also make 
precision measurements at hadron collider experiments, collab¬ 
oration between PDF-generators and experimenters is essen¬ 
tial. The W mass illustrates where this would be useful prob¬ 
ably better than any other precision electroweak measurement. 
Given that the future seems to be evolving towards higher en¬ 
ergy hadron colliders, the necessity of errors on parton distri¬ 
bution functions can only increase, as will the care required in 
using these functions correctly. 


V. HEAVY QUARK HADROPRODUCTION 


Improved experimental measurements of heavy quark ha- 
droproduction have increased the demand on the theoretical 
community for more precise predictions. [^1]] The first Next-to- 
Leading-Order (NLO) calculations of charm and bottom hadro- 
production cross sections were performed some years ago.|32| 
As the accuracy of the data increased, the theoretical predictions 
displayed some shortcomings: 1) the theoretical cross-sections 
fell well short of the measured values, and 2) they displayed 
a strong dependence on the unphysical renormalization scale 
/r. Both these difficulties indicated that these predictions were 
missing important physics. 

One possible solution for these deficiencies was to consider 
contributions from large logarithms associated with the new 
quark mass scale, such as^ ln(s/mQ) and ln(py/TOQ), Push¬ 
ing the calculation to one more order, formidable as it is, would 
not improve the situation since these large logarithms persist to 
every order of perturbation theory. Therefore, a new approach 
was required to include these logs. 

In 1994, Cacciari and Greco [p4|] observed that since the heavy 
quark mass played a limited dynamical role in the high pt re¬ 
gion, one could instead use the massless NLO jet calculation 
convoluted with a fragmentation into a massive heavy quark pair 
to more accurately compute the production cross section in the 
region pt mg. In particular, they find that the dependence on 
the renormalization scale is significantly reduced, (cf.. Fig. p^ . 

A recent study [|3^ investigated using initial-state heavy quark 
PDFs and final-state fragmentation functions to resum the 
large logarithms of the quark mass. The principle ingredi¬ 
ent was to include the leading-order flavor-excitation (LO- 
FE) graph (Fig. and the leading-order flavor-fragmentation 
(LO-FF) graph (Fig. E> in the traditional NLO heavy quark 
calculation, These contributions can not be added naively to 


°Here, ttiq is the heavy quark mass, s is the energy squared, and pjp is the 
transverse momentum. 



pp ^BX, Vs =1.8TeV 



Figure 16: Heavy quark hadroproduction data. Cf., Ref. 
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Figure 17: Scale dependence of the heavy quark hadropro- 
duction cross section as a function of p. = ^Pre/ at y = 0 
and pt = 80 GeV. The NDE curve is the calculation of 
Ref. The fragm., funct. and bom curves are the calcu¬ 

lation of Ref. [Q. 



Figure 18: a) Generic leading-order diagram for flavor- 

excitation (LO-FE), gQ gQ. b) Subtraction diagram for 
flavor-excitation (SUB-EE), ® <y{gQ gQ)- c) Next- 

to-leading-order diagram for flavor-creation (NLO-EC). 



Eigure 20: The differential cross section (Pa/dp'^/dyi at 
PT = 20, 40GeU and yi = 0 in {pb/GeV"^) vs. p. The lower 
curves (thin line) are the heavy quark production cross sections 
ignoring flavor-excitation (EE) and flavor-fragmentation (EE). 
The upper curves (thick line) are the heavy quark production 
cross sections including FE and EE. Cf, Ref. [^]. 


the 0{a^g) calculation as they would double-count contributions 
already included in the NLO terms; therefore, a subtraction term 
must be included to eliminate the region of phase space where 
these two contributions overlap. This subtraction term plays the 
dual role of eliminating the large unphysical collinear logs in the 
high energy region, and minimizing the renormalization scale 
dependence in the threshold region. The complete calculation 
including the contribution of the heavy quark PDEs and frag¬ 
mentation functions 1) increases the theoretical prediction, thus 
moving it closer to the experimental data, and 2) reduces the 
p-dependence of the full calculation, thus improving the pre¬ 
dictive power of the theory. (Cf., Eig p0[ ) 

In summary, heavy quark hadroproduction is of interest ex¬ 
perimentally because of the wealth of data allows precise tests 
of many different aspects of the theory, namely radiative correc¬ 
tions, resummation of logs, and multi-scale problems. Hence, 
this is a natural testing ground for QCD, and will allow us to ex¬ 
tend the region of validity for the heavy quark calculation. This 
is an essential step necessary to bring theory in agreement with 
experiment. 


VI. SUMMARY 

A. Kinematic Reach of Future Machines 


, , Table V: Euture ep collider machines chosen for study. 


Eigure 19: a) Generic leading-order diagram for flavor- 

fragmentation (LO-EE), a{gg gg) G Dg^q. b) Subtraction 
diagram for flavor-fragmentation (SUB-EE), cr{gg gg) 0 
^dg^Q. c) Next-to-leading-order diagram for flavor-creation 
(NLO-EC). 


Index 

Elepton 

(Gev) 

Eproton 

(Gev) 

(Gev) 

Machine(s) 

1 

27 

820 

300 

Hera 

2 

35 

7,000 

990 

Lep X LHC 

3 

8 

30,000 

980 

Low E lepton x 60 GeV pp 

4 

30 

30,000 

1900 

Lep X 60 GeV pp 

5 

500 

500 

1000 

NLC X conv. p 

6 

2,000 

500 

2000 

p collider x conv. p 
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A central goal of this workshop was to study the physics po¬ 
tential of future facilities. Here, we focus on lepton-hadron col- 


































Figure 21: Kinematic reach of future facilities. 



Figure 22: Kinematic reach of future facilities. 

liders. We expand our study beyond the single ep machine pro¬ 
posed in the workshop outline, and consider a mix of lepton 
and hadron beams from those proposed for the lepton-lepton 
and hadron-hadron options. The complete list is given in Ta¬ 
ble ^ To covert these parameters into the {x, Q^} range, we 
make use of: 


2/ = 1 - ^(1-cos6»f) , 

(7) 

= 2EeE'^{l + cos Bg) , 

(8) 

and 

x = Q^lsy . 

(9) 

For collider kinematics, we use 

S ~ AEeEp . 

(10) 


Here, is the incoming lepton energy, E'^ is the outgoing lep¬ 
ton energy, Ep is the incoming hadron energy, and 6^ is the 
lepton scattering angle. 

To set practical limits on measurement of the final state, we 
impose: 


LHC 



Figure 23: Kinematic reach of present and planned facilities. 
Note the full {a:, Q^} region is clipped by the plot. 

• y > 0.01 (resolution), 

• y < 1 (kinematic cut) 

• 9e> 10° 

• Be < 179° 

The constraint B^ < 179° may be somewhat optimistic; if we 
relax this to 9g ;$ 176°, the result is to lose some of the low Q 
region. The constraint Bg > 10° has a relatively small effect; 
for the higher energy machines {e.g, 2 & 3), it clips the upper Q 
region. 

We display the kinematic reach for these proposed machines 
in Figs. ^ and ^ We include HERA for reference. In Fig. ^ 
we show the three machine options with a CMS of ^/s ~ 1 TeV. 
In Fig. ^ we show HERA and the remaining two machine op¬ 
tions. In Eig. ^ we show the present and planned (LHC) facil¬ 
ities. 

Although there is currently no plan to extract the primary 
beam to make a neutrino fixed target experiment at either the 
LHC or at a 2 TeV muon collider, there is a case to be made for 
doing precisely that. First of all, it would be very interesting to 
see if there were an anomalous rise in xF^ similar to that seen in 
F2 at HERA. Secondly, the low x region of the Bjorken integral 
is anomalously large, and an outstanding question is, what is 
the very low x behavior of the Gross-Llewellyn Smith integral 
{xF^)l Either an experiment at the LHC or one at a 2 TeV muon 
collider could extend the range of the ’’Eixed Target” region in¬ 
dicated in Eigure 23 by an order of magnitude in the log (1/a;) 
direction, assuming an order of magnitude higher neutrino en¬ 
ergies than what CCER/NuTeV has. The neutrino cross section 
would be an order of magnitude higher than the one applicable 
for CCER/NuTeV, so good statistics are in principle attainable. 
Although these experiments would not have the kinematic reach 
to extremely low x that ep machines have, they can measure 
to high precision the non-singlet structure function, which at 
present has only been measured down to a; = .01. In principle 
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an ep machine running with both positive and negative leptons 
could do the same, but the luminosity requirements may be pro¬ 
hibitively high. We have still not learned all that we can learn 
from neutrino experiments, and even modest improvements in 
neutrino energies can uncover much new ground. 

While we would of course like to probe the full {x,Q^} 
space, there are some particular reasons why the small x re¬ 
gion is of special interest. For example, the rapid rise of the 
F 2 structure function observed at HERA suggests that we may 
reach the parton density saturation region more quickly than an¬ 
ticipated. Additionally, the small x region can serve as a useful 
testing ground for BFKL, diffractive phenomena, and similar 
processes. We can clearly see in Fig. E that with a hxed 
we can best probe the small x region with a high energy hadron 
beam colliding with a low energy lepton beam, and the loss in 
the high Q region is minimal. From these (preliminary) studies, 
it would seem the optimal ep facility would match the high¬ 
est energy hadron beam available with a modest energy lepton 
beam. 
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